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on one-loop coefficients with that of mean-field improvement of the operators. 
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1. Introduction 

Matrix elements of electroweak four-fermion operators between hadronic states play a central 
role in constraining the parameters of the CKM matrix. The Wilson coefficients of the operators 
are typically calculated in a continuum renormalization scheme (usually the NDR [Naive Dimen- 
sional Regularization] scheme with MS subtraction), so one needs to convert the matrix elements 
calculated on the lattice into this continuum scheme. The multiplicative factors that are needed 
(which are, in general, matrices) are called matching factors or "Z-factors". We consider here the 
matching factors for four-fermion operators constructed using HYP-smeared staggered fermions 
[jl|], and calculate them at one-loop order in perturbation theory. These are needed for our ongoing 
calculation of Bk using HYP-smeared staggered fermions [§] on the MILC "asqtad" lattices 

The new feature compared to the previous calculation [Q] of four-fermion matching factors 
is the use of the Symanzik-improved gluon action instead of the Wilson plaquette action. This 
extension is necessary because the MILC lattices are generated using the former action. 

2. Action and Feynman Rules 

For valence quarks we use the HYP-smeared action, which is simply the unimproved staggered 
action 

L n,n n 

(where Tfy (n) = (— j an d we have set the lattice spacing a to unity) except that the links, 
, are HYP-smeared fat links. By contrast, in the unimproved action the links are replaced by the 
original, "thin", links, U^. The relation between these two types of links is given in Ref. [jl|]. What 
we need is the relation in perturbation theory and in momentum space: 

V M (n) = exp(igfl M (n)) (2.2) 
V») = /^(04l^vWAvW^ ( ' ,+/i/2) + ^(A 2 ). (2.3) 

Here A v are the momentum space gluon fields obtained from the unimproved links, and h llv {k) is 
a smearing kernel containing all relevant information about the HYP smearing. For the HYP links 
whose coefficients are chosen to remove &(dr) taste symmetry breaking at tree level (a\ = 0.875, 
a 2 = 4/7 and a 3 = 0.25 in the notation of Ref. [|I|]), h ilv (k) is [g, g, [7|] 

h^iv(k) = 8^ v D^(k) + (1 - d^s^SvGv^ik) , (2.4) 

where 

D M (*) 
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Here, ^/v/p/a and = sm(k^/2). We recover the unimproved action (in which the links 
are thin) by setting h iiV {k) = 8^ v . 

One would expect the @{A 2 ) term in Eq. ( j2j| ) to be needed at one-loop order, since it can 
give rise to tadpole diagrams. It turns out, however, that its contribution vanishes due to the SU (3) 
projection used in HYP-smearing, as discussed in Refs. ||, g]. 

We use the tree level Symanzik-improved gluon action [|l(], which can be written as 



S8 ~8l 



|£ReTr(l-f/ pl ) - i-£ReTr(l - U«) 

-'pi 1Z rt 



(2.7) 



where "pi" and "rt" represent plaquette and rectangle, respectively. In fact, the MILC collaboration 
use a (partially) one-loop improved Symanzik gluon action to generate their configurations [|l2|], 
and not the tree-level form. However, the one-loop corrections in the gluon action do not impact 
the matching factors of fermionic operators until two-loop level, and thus it is sufficient to use 
the tree-level gluon action. Similarly, the effect of dynamical quarks does not enter until two-loop 
level, so the fact that the valence and sea quarks have different actions does not impact our one-loop 
calculation. 

Feynman rules for the staggered action and four-fermion operators can be found in literature, 



and we do not repeat them here. Quark propagators and vertices can be found in Refs. [13, 14, 
|8p. The generalization to various improved staggered fermions such as asqtad and HYP-smeared 
staggered fermions can be found in Ref. |]]. The improved gluon propagator can be found in 
Refs. 



3. Four-fermion Operators 



We use the hypercube construction of Ref. [ ]15[ ] for our staggered four-fermion operators. 
These come in two kinds, which differ in the contractions of their color indices. First we have 
"one color-trace" operators, labeled with a subscript /, whose form is 

[5xF][5'x4(y) = i I lX { a l \2y+A)J^r) AB% l 2) (2y + B)] 

4 A,B,A',B' 

x [£ ] (2y +A')(yy^W*,? ) + #)] 
xy ah '(2y+A,2y + B')y a ' b (2y+A',2y + B). (3.1) 

Here y £ Z 4 is the coordinate of hypercubes, hypercube vectors S and S' denote the spins of the 
component bilinears, while F and F' denote the tastes. Indices a, b, a', and b' denote colors, while 
superscripts (i) for i = 1,2,3,4 label different flavors (not tastes). Different flavors are chosen to 
forbid penguin diagrams. Two fat Wilson lines T ab ' (2y +A,2y +B') and f a ' b (2y +A',2y + B) 
ensure the gauge invariance of the four-fermion operators. A fat Wilson line ^ /afo (2y + A, 2y+B'), 
for example, is constructed by averaging over all the shortest paths connecting 2y +A and 2y + B', 
with each path formed by products of HYP-smeared links V^. When we use the unimproved 
staggered action the Wilson lines are composed of unsmeared links, U^. 
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The second kind of operator is called "two color-trace", and takes the form 



[SxF][S'xF'] II (y) = i- 4 I [xi l \2y+A)(r s 0^) AB xl 2> (2y + B)} 



,( 2 ), 



44 



A,BA',B' 



x \x£ ] (2y +A')( Ys >® Zf>) aV xP (2y + B')] 
x Y ab (2y+A,2y + B) y a ' b ' (2y+A', 2y + B') 



(3.2) 



where the subscript // denotes two color-traces. These operators differs from those with one color- 
trace by the choice of fat Wilson lines — here they connect within each bilinear, whereas for the one 
color-trace operators they connect between bilinears. 

Following Refs. [||, [14], [l^], we also consider mean-field improved operators. For the unim- 
proved action and operators this is usually referred to as tadpole improvement [17]. For the HYP- 
smeared action and operators, mean-field improvement is achieved by rescaling the staggered fields 
and the link matrices: 



X^Y 



voX, 



X^Y 



voX, 



Vii =V^/v , 



(3.3) 



where vq is the mean-field improvement factor defined by 



v = 



r ReTr(y pl ) 



1/4 



(3-4) 



with Vpi the plaquette composed of HYP-smeared links. 



4. Renormalization of Four-fermion operators 

Analytic formulae for one-loop perturbative corrections to the four-fermion operators for 
HYP-smeared staggered fermions with the Wilson plaquette action are given in Ref. Jjjj. It turns 
out that the generalization to using the improved gluon action rather than the Wilson plaquette 
action can be achieved simply by replacing the "composite" gluon propagator; 

{\/k 2 )Y J h^h vX ^Y. h ^h v ^ a p, (4-1) 

where 3> a p is the improved gluon propagator, while the Wilson gluon propagator is (8 a p/k 2 ) with 
k 2 = 4£ Q ,sin 2 (& a /2). This simplification holds both for bilinear operators [^] and four-fermion 
operators. 

The Feynman diagrams which contribute to one-loop matching factors for four-fermion opera- 
tors can be found in Ref. [H]. We have undertaken two independent calculations, using significantly 
different methods, to cross-check our results. Another check is obtained by using Fierz transfor- 
mations, which interchange one and two color- trace operators. The form of these transformations 
relevant to our set-up is 

(ys®&) A b{7s>®Zfi)a>b> = 77 L {ysyl®^F<) AB ,{ys>yM®^F) A , B , (4.2) 

10 MN 

where M, N, A, A', B, and B' are hypercube vectors. 
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The result of the calculation can be expressed as 

d Cont,(l) (At) = ^ z .. (At)fl)5 Lat (1/a)) ^ 

j 

where the superscripts "Cont" and "Lat" represent the continuum and lattice respectively. The sub- 
indices i and j run over all combinations of spins and tastes of the four-fermion operators that are 
allowed by the lattice symmetries. The vectors on the four-fermion operators indicate that one and 
two color-trace operators, for given spins and tastes, are collected into a two-dimensional vector. 
We call the vector space in which these vectors live the "color-trace" space. The renormalization 
scale in the continuum scheme is }X, and a (now reintroduced) is the lattice spacing. Finally, the 
superscript "(1)" indicates that the continuum operator has been matched at one-loop order. 
The matching factor Z,- ; takes the form 



'J ~ "'J ^ (471)2 



fijlog(na) + 



Ci 



(4.4) 



where fjj and en are, respectively, the anomalous dimension matrix and the finite coefficients. 
Apart from their explicit indices, both are also matrices in the color-trace space, as denoted by the 
"hats". The finite coefficients are given by the difference of finite terms in the continuum and lattice 
one-loop calculations, 



cf? Dt -c}f. (4.5) 



5. Matching Factors for Bk 

We have obtained the matching factors for all four-fermion operators of the forms given above, 



i.e. for all choices of S, S', F and F'. Detailed results will be presented elsewhere [ 18 1. We consider 
here only the matching factors needed for Bk- The continuum four-fermion operator relevant to Bk 
is 

&b? = [^7m(1 " 7 5 K'][^(1 - Ys)d b ] . (5.1) 



£^ = y I + y [I + ^ [ + £ / in (5.2) 



The corresponding operator on the lattice, in the "two spin-trace" formulation Q19j |20| ], can be 
written at tree-level as 

4 

yhat \ /?Lat 

i=l 

where 

^Lat = y [= ^ x p] ^ x p][ ; ^>Lat ^ % = ^ x p] ^ x p][[ ^ (5.3) 

^ 3 Lat = srfj = [A M x P] [Ap xP] r , =^ II = [A tl x P] [A„ xP] n . (5 .4) 

Here, V^, A^, and P stand for the vector, axial- vector, and pseudo-scalar respectively. The taste of 
all the operators are taken to be the pseudo-scalar P, i.e. F = ^5. 

At one-loop order, many lattice operators contribute to the matching formula. In practice, 
however, in our numerical calculation of Bk we keep only those operators in which the bilinears 
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have taste It turns out that only the four operators listed above have this taste. Thus what is 
relevant for our numerical calculation is the truncated one-loop matching formula 



j)Cont, (1), Trunc 



s-Lat, (1) 



1 + 



(An) 



- A log {lxa) + \ 



11 

T 



Lat 



(5.5) 



— for all four 



We use here the NDR scheme in the continuum, for which the finite coefficient is 
operators. 

Results for c\ ai are given in table H] for various choices of lattice operator. We see that, in most 
cases, the coefficients are slightly reduced by improving the gauge action, although the reduction 
is smaller than that obtained by for HYP fermions by mean-field improvement. As can be seen 
from Eq. (5.5), however, reducing the c^ at does not necessarily lead to matching factors closer to 



unity. A better measure is the range of the matching coefficients, i.e. the largest difference between 
them, for the continuum and logarithmic contributions then cancel. This quantity is also given in 
the Table, and one sees that it is moderately reduced by improving the gauge action. 





(a) 


(b) 


(c) 


(d) 


to 


(/) 


Gluon action 


Wilson 


Sym 


Wilson 


Sym 


Wilson 


Sym 


Quark action 


Naive 


Naive 


HYP 


HYP 


HYP 


HYP 


Mean-field impr.? 


Y 


Y 


N 


N 


Y 


Y 


„Lat 
c l 


-2.35 


-2.49 


-4.98 


-3.65 


-2.17 


-1.72 


„Lat 
°2 


-12.92 


-11.53 


-11.17 


-8.58 


-5.49 


-4.73 


Xat 
6 3 


-2.95 


-3.08 


-5.50 


-4.12 


-2.69 


-2.19 


Xat 

c 4 


-3.73 


-2.90 


1.01 


1.09 


1.01 


1.09 


Range 


10.57 


9.04 


12.18 


9.67 


6.50 


5.82 



Table 1: Results for c^ 31 for various choices of gauge and fermion action. We quote only two decimal places 
for brevity. 

To give a sense of the size of the matching coefficients themselves, we show in table || re- 
sults for the ultrafine MILC lattices (a « 0.045 fm), setting jl = \ /a ("horizontal matching"), and 
g 2 /(An) = 0.2096 (the value in the MS scheme at pL = \/a). For the actions we use in practice 
[column (d)] the one-loop corrections range between ±8%. 





(a) 


(b) 


to 


(d) 


to 


(/) 




0.978 


0.980 


1.022 


1.000 


0.975 


0.968 


C2 


1.154 


1.131 


1.125 


1.082 


1.030 


1.018 


C3 


0.988 


0.990 


1.031 


1.008 


0.984 


0.975 


c 4 


1.001 


0.987 


0.922 


0.921 


0.922 


0.921 



Table 2: Values of c, for the MILC ultrafine lattice. Notation for actions as in table U 



'The rationale for this is that we use external kaons with taste §5. As shown in Ref. pl|], however, leaving out 
the operators with other tastes leads to an error of ^(c^m^/Ay), which is of next-to-leading order in staggered chiral 
perturbation theory. This error must be accounted for when fitting. 
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6. Conclusion 

We have calculated the matching factors for four-fermion operators using various fermion 
and gauge actions. Most useful are our results for the fermion action and operators containing 
HYP-smeared links, and with the gluon action being Symanzik improved, for these are needed 
for our parallel calculation of Bk, and have been used in Ref. For these choices, the one- 
loop corrections are of moderate size, with the range of corrections being m 10 x a s /(4n) ph a s , 
which is the naively expected size. The impact of improving the gauge action turns out to be 
numerically small, much less than the reduction in the size of the corrections caused by using 
mean-field improved operators. 
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